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ABSTRACT 

We have made velocity observations of the star a Cen B from two sites, allowing us to identify 37 oscillation 
modes with I = 0-3. Fitting to these modes gives the large and small frequency separations as a function of 
frequency. The mode lifetime, as measured from the scatter of the oscillation frequencies about a smooth trend, 
is similar to that in the Sun. Limited observations of the star d Pav show oscillations centred at 2.3 mHz with 
peak amplitudes close to solar. We introduce a new method of measuring oscillation amplitudes from heavily- 
smoothed power density spectra, from which we estimated amplitudes for a Cen A and B, /3 Hyi, S Pav and the 
Sun. We point out that the oscillation amplitudes may depend on which spectral lines are used for the velocity 
measurements. 

Subject headings: stars: individual (a Cen A, a Cen B, (3 Hyi, 6 Pav) — stars: oscillations — 
Sun: helioseismology 



1. INTRODUCTION 

The a Cen A & B system is an excellent target for asteroseis- 
mology. Velocity oscillations in the A component were mea- 
sured by Bouchy & Carrier (2001, 2002) and subsequently by 
Butler et al. (2004, see also Bedding et al. 2004). The first de- 
tection of oscillations in a Cen B (HR 5460) was made by Car- 
rier & Bourban (2003), based on velocity measurements with 
the CORALIE spectrograph in Chile spanning 13 nights. They 
observed excess power centred at 4 mHz and used the autocor- 
relation of the power spectrum to measure a large separation of 
161.1 yuHz. They identified twelve oscillation frequencies with 
I = 0-2, although the high sidelobes from their single-site ob- 
servations, coupled with the relatively low signal-to-noise (2.5- 
3.5), means that these may have been affected by daily aliases. 
Theoretical modelling of a Cen A and B has been carried out 
several times, most recently by Eggenberger et al. (2004), who 
also give a thorough review of previous work. 

Here we report observations of a Cen B made from two sites 
that show oscillations clearly and allow us to measure their fre- 
quencies, amplitudes and mode lifetimes. 

2. VELOCITY observations AND POWER SPECTRA 

We observed a Cen B in May 2003. At the European South- 
ern Observatory in Chile we used UVES (UV- Visual Echelle 
Spectrograph) at the 8.2-m Unit Telescope 2 of the Very Large 
Telescope (VLT)''. At Siding Spring Observatory in Australia 
we used UCLES (University College London Echelle Spectro- 
graph) at the 3.9-m Anglo-Australian Telescope (AAT). In both 
cases, an iodine cell was used to provide a stable wavelength 
reference (Butler et al. 1996). 

At the VLT we obtained 3379 spectra of a Cen B, with typi- 
cal exposure times of 4 s and a median cadence of one exposure 
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every 32 s. At the AAT we obtained 1642 spectra, with typi- 
cal exposure times of 10-12 s (but sometimes as long as 30 s in 
poor weather) and a median cadence of one exposure every 91s. 
Note that, unlike for our observations of a Cen A, UCLES was 
used in standard planet-search mode. For a Cen A we rotated 
the CCD by 90 degrees to speed up readout time, but found that 
this introduced drifts and sudden jumps in the velocities that 
we believe to be related to the movement of the CCD dewar as 
liquid nitrogen boiled off and was refilled (Butler et al. 2004). 

The resulting velocities are shown in Fig. 1, and the effects 
of bad weather can be seen (we were allocated five nights with 
the VLT and eight with the AAT). Both sets of velocities show 
trends during each night that are presumably due to a combi- 
nation of instrumental drift and stellar convection noise. To 
remove these trends, we have subtracted a smoothed version of 
the data one night at a time, producing the detrended time se- 
ries shown in Fig. 1. This process of high-pass filtering is not 
strictly necessary, but it does slightly reduce the noise leaking 
into higher frequencies that would otherwise degrade the oscil- 
lation spectrum. 

Note that Fourier analysis on unevenly spaced data cannot 
use the Fast Fourier Transform (EFT) algorithm. We instead 
employed the method used for many years by various groups, 
including our own: we calculated the discrete Fourier trans- 
form, but with individual data points being weighted according 
to their quality (e.g. Deeming 1975; Frandsen et al. 1995). 

For about one hour at the end of each night, when a Cen B 
was inaccessible, we observed the star (5 Pav (HR7665; 
HIP99240; G8 V). We obtained a total of 179 spectra with the 
VLT (median cadence 58 s) and 77 with the AAT (median ca- 
dence 181s). The velocities are shown in Fig. 2 and there is a 
clear periodicity of about 7 min. We did not expect to be able 
to measure frequencies with such a limited data set, but we do 
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see a clear excess in the power spectrum from the oscillations 
(Fig. 3), whose amplitudes we discuss in Sec. 5. 

Our analysis of the velocities for a Cen B follows the method 
that we developed for a Cen A (Butler et al. 2004). We have 
used the measurement uncertainties, ct,, as weights in calculat- 
ing the power spectrum (according to w, = 1 jaf), but modified 
some of the weights to account for a small fraction of bad data 
points. In this case, only three data points from UVES and 
none from UCLES needed to be down- weighted. After these 
adjustments, we measured the average noise in the amplitude 
spectrum at high frequencies (above the stellar signal) to be 
1.49cms-i for UVES (6-8 mHz) and 4.28 cm s"! for UCLES 
(4.75-5.50mHz). Note that we used a lower frequency range 
to measure the noise in the UCLES data because the sampUng 
time of those observations was significantly longer than for 
UVES (see above). 

The power spectrum of the combined series is shown in the 
top panel of Fig. 4. We refer to this as the noise-optimized 
power spectrum because the weights have been chosen to min- 
imize the noise. The noise level in amplitude is 1.39 cm s"' (6- 
8 mHz). As such, these measurements replace our observations 
of a Cen A (L9 cm s~^) as the most precise amphtude spectrum 
obtained on any star apart from the Sun. 

The middle panel of Fig. 4 shows a close-up of the noise- 
optimized power spectrum. There is a clear series of regular 
peaks with a spacing of about 81 /xHz, and we therefore con- 
firm that the large separation is about 162 /iHz, in agreement 
with the value reported by Carrier & Bourban (2003). The inset 
shows the spectral window, in which we see prominent side- 
lobes (38% in power) due to gaps in the observing window. 

As for a Cen A, we have also generated a power spectrum in 
which the weights were adjusted on a night-by-night basis in or- 
der to to minimize the sidelobes. The result, which we refer to 
as the sidelobe-optimized power spectrum, is shown in the bot- 
tom panel of the figure. The values for these weight multipUers 
were: 1.0 for all UCLES nights except the last, which had 0.8; 
and 3.2, 3.2, 2.3 and 2.3 for the four UVES nights. Adjusting 
the weights in this way increased the noise level to 2.40 cms"' 
(6-8 mHz), but allowed us to identify correctly peaks that might 
otherwise be obscured by sidelobes from their stronger neigh- 
bours (now reduced to 13% in power). In addition, the side- 
lobe minimization has slightly improved the frequency reso- 
lution, with the FWHM of the spectral window decreasing by 
about 20% (from 1.83 /xHz to 1 .44 /xHz). This is because of the 
increased weight given to the UCLES data, which covers the 
longer timespan. Note that the dotted lines in Fig. 4 show our 
final oscillation frequencies, which are discussed in detail in the 
next section. 

3. OSCILLATION FREQUENCIES 

We measured the frequencies of the strongest peaks in the 
power spectrum in the standard way, using iterative sine-wave 
fitting. We did this for both the noise-optimized and the 
sidelobe-optimized power spectra, and the resulting frequen- 
cies are plotted in echelle format in Fig. 5. The echelle for- 
mat takes advantage of the fact that mode frequencies for low- 
degree p-mode oscillations are approximated reasonably well 
by the asymptotic relation: 

= Az/(n + 5/ + e)-/(/+l)Oo. (1) 
Here n (the radial order) and I (the angular degree) are integers, 
IS.V (the large separation) reflects the average stellar density, ZJq 
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is sensitive to the sound speed near the core and e is sensitive to 
the surface layers. 

Symbol size in Fig. 5 indicates the signal-to-noise ratio 
(SNR) of the peaks, and all peaks with SNR > 2.5 are shown. 
Many of these represent oscillations but some are artefacts due 
to the non-linear nature of the iterative fitting method and its in- 
teraction with noise. The dotted lines represent a fit to the final 
frequencies that is discussed below; in this figure, they serve 
to guide the eye and allow us to identify ridges with angular 
degrees of Z = 0, 1,2 and 3. 

We can see from Fig. 5 that many peaks were identified in 
both versions of the power spectrum and, not surprisingly, many 
were only above SNR = 2.5 in the noise-optimized version. 
However, four peaks that were only significant in the sidelobe- 
optimized version lie close to the oscillation ridges (three with 
/ = and one with / = 1). The detection of these peaks vin- 
dicates our decision to examine the sidelobe-optimized power 
spectrum. 

The next step was to select those peaks that we believe are 
due to oscillations and reject those due to noise. This is the 
only subjective part of the process, but is required if we are 
to measure as many frequencies as possible from our data be- 
cause at this SNR, not all the extracted peaks will be genuine. 
Our final set of identified modes is shown in Fig. 6 and Table 1 
(and also indicated by the dotted lines in Fig. 4). In some cases 
the iterative sine-wave fitting produced two peaks (in one case, 
three) that appear to arise from a single oscillation mode, which 
is to be expected if we have partially resolved the modes (see 
Sec. 4 for discussion of the mode hfetimes). In these cases, 
and in cases where the same peak was detected in both versions 
of the power spectrum, all peaks are shown in Fig. 5 but they 
are combined into a single unweighted mean in Fig. 6 and Ta- 
ble 1 (the uncertainties given in the table are discussed below in 
Sec. 4). We have been conservative in not selecting three peaks 
outside the main region because we cannot be sure of the cur- 
vature of the extrapolated ridge lines. These are: 2827.7 /uHz 
(/ = 0?), 5894.9 /iHz (/ = 2?) and 5973.1 /xHz {I = 1?). 

Inspection of Fig. 4 shows some apparent mismatches be- 
tween the peaks in the power spectrum and the dotted Unes rep- 
resenting the extracted frequencies. For example, there are en- 
hancements in power at 3600, 3700 and 4050 /iHz that do not 
correspond to identified modes. We noted similar structures in 
our analysis of a Cen A (Bedding et al. 2004) and the explana- 
tion is the same: the interaction of the window function with 
noise peaks and oscillation peaks. We have verified that the 
same phenomenon occurs in solar data by analysing segments 
of the pubhcally-available^ 805-day series of full-disk veloc- 
ity observations taken by the GOLF instrument (Global Oscil- 
lations at Low Frequencies) on the SoHO spacecraft (Ukich 
et al. 2000), which have a samphng time of 80 s. We applied 
the same sampling window as our observations of a Cen B and 
the results showed low-level structure in the power spectrum 
similar to those in Fig. 4. This confirms that such a structure is 
a natural consequence of the spectral window interacting with 
multi-mode oscillations having finite lifetimes. 

In Fig. 6 the open squares show the frequencies reported by 
Carrier & Bourban (2003). It is clear that, while they found the 
correct large separation, there is a shift in their frequencies of 
one cycle per day (11. 57 /iHz). This reflects the problem asso- 
ciated with single-site observations, especially when the peaks 
have low signal-to-noise. 
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As with a Cen A, we note a scatter of the peaks about the 
ridge lines that is much higher than expected from SNR con- 
siderations and which we interpret as being due to the finite 
lifetime of the modes. We therefore fit to the ridges, in order 
to obtain more accurate estimates for the eigenfrequencies of 
the star. In the case of a Cen A (Bedding et al. 2004) we fit- 
ted to the frequencies in a two-step process by first fitting the 
three small separations and then fitting a parabola to the indi- 
vidual frequencies. Here, we adopted a simpler approach that 
gives almost identical results, in which we fitted directly to the 
frequencies (see Ulrich 1986). The nine fitted parameters spec- 
ify: the curvatures of the parabolas (one common value), the 
large separation at some reference frequency for each / (four 
values) and the absolute positions of the ridges (four values). 
The equations for this fit are as follows: 

v„f^= [3950.57 + 161. 45 (n - 23) + 0. 101 (n- 23)^] /iHz (2) 

[4026.23+ 161.28 (n-23) + 0.101(«-23)2] /iHz (3) 

i/„2= [4101.41 + 161.63 («-23) + 0.101(n-23)^] [iVlt (4) 

jy„3= [4171.13+ 161.76(n-23) + 0.101(n-23)^] ^Hz. (5) 

We show this fit as the dotted curves in Figs. 5 and 6. 

Figure 7 shows the small frequency separations (top panel), 
the £)o parameter (middle panel) and the large separations (Ai/; 
bottom panel), using the same definitions as Bedding et al. 
(2004). Thus, (5t'()2 is the separation between adjacent peaks 
with I = and I = 2 and Si^u is the separation between I = 1 
and 1 = 3. The third small separation, Svqi, is the amount by 
which 1 = 1 modes are offset from the midpoint between the 
/ = modes on either side: 

5uoi = j {Unfl + iy„+ifl) -t^n,i- (6) 
Since one could equally well define Si^m to be the offset of 1 = 
modes from the midpoint between consecutive 1=1 modes, 

Sm = '^n,0-j{'^n-l,l+'^n,l), O) 

we have shown both versions in Fig. 7 (upper panel). Note that 
Do, which is sensitive to the stellar core, corresponds to ^t^foi, 

51^13 and ^Svoi and is a constant if the asymptotic relation 
holds exactly. The solid Unes in the top and bottom panels show 
the separations calculated from the fitted equations (in the bot- 
tom panel, only Aj/q is shown; lines for other I values are almost 
indistinguishable). The upward trend of Av with frequency is 
responsible for the curvature in the echelle diagram. 

In Table 2 we give the large and small frequency separa- 
tions for a Cen B at a frequency of 4.0 mHz. We also give e, 
which is a dimensionless quantity commonly used to param- 
eterize the absolute position of the frequency spectrum - see 
equation (1). We are not able to detect any statistically signif- 
icant variations in the large separation with angular degree (1); 
the weighted mean value ( A;^) has an uncertainty of 0.04%. We 
also find that the value of Do (as defined above) at 4.0 mHz is 
the same within uncertainties for all three small separations, 
and our mean value for Do has an uncertainty of 3%. We 
can therefore place a Cen B in the so-called asteroseismic H- 
R diagram (Christensen-Dalsgaard 1984; see also Ulrich 1986; 
Gough 2003; Oti Floranes et al. 2005), and these frequency sep- 
arations should be compared with theoretical models. 

Detailed modelling of the a Cen system was carried out by 
Eggenberger et al. (2004). They determined the parameters of 
the system through a least-squares fit to the observed quantities, 
including both 'classical' photomettic and spectroscopic data 
and oscillation frequencies from Bouchy & Carrier (2002) and 



Carrier & Bourban (2003). Their two preferred models (Mod- 
els Ml and M2) both have an age of 6.5 Gyr. For a Cen B they 
obtained large separations A^o of 161.7 and 161.1 /xHz, respec- 
tively, while the / = to 2 small separations were ^^"02 = 10.3 
and 10.2 /iHz, respectively; both are essentially consistent with 
the values found here (cf. Table 2), although with some slight 
preference for Model Ml. It should be noted, however, that the 
average values of these quantities depend on the detailed way 
in which the averages are computed, including the selection of 
modes. Also, the computed values of At'o, and to a lesser extent 
6i'o2, are affected by the uncertain physics of the near-surface 
layers of the model. An analysis of the combined set of fre- 
quencies for a Cen A, presented by Bedding et al. (2004), and 
the results of a Cen B shown in Table 1 is in progress. 

3.1. Curvature at high and low frequencies 

We now describe a method that allows us to measure the os- 
cillations at low SNR, outside the central region in which we 
can identify individual modes. This allows us to examine cur- 
vature in the echelle diagram, which corresponds to measuring 
the large separation as a function of frequency. The method 
relies on smoothing the power spectrum in order to increase 
the contrast between the oscillation signal and the background 
noise. The first step was to smooth the power spectrum quite 
heavily, with a FWHM of Az//4. This type of smoothing was 
used by Garcia et al. (1998) to find high-frequency peaks in the 
solar power spectrum, but here we use an even broader smooth- 
ing function. The smoothing combines the power from each 
even-degree pair of modes (/ = and 2) into a single resolution 
element, and the same applies to the odd-degree pairs (/ = 1 
and 3). 

We next arranged the smoothed power spectrum in echelle 
format, with frequencies reduced modulo 5 Aj/ rather than the 
conventional Ai/. This causes power associated with odd and 
even pairs to line up, allowing us to smooth in the vertical direc- 
tion in this echelle diagram to improve the contrast further (for 
this smoothing we used FWHM = 3.4Az/). Finally, we mea- 
sured the highest peak in each half-order and these are marked 
by the open circles in Fig. 8. We see a clear ridge of power 
that extends well beyond the central region. The dotted lines 
are the fits shown in Fig. 6 and given by equations (2)-(5). We 
see good agreement in the central region, but the smoothed data 
allow us to measure curvature well outside the region in which 
we were able to identified individual oscillation modes. 

In order to evaluate this method, we have also applied it to 
the oscillations in the Sun by analysing the 805 -day GOLF se- 
ries (see Sec. 3). We used the same amount of smoothing as for 
a Cen B (when specified in terms of Au) and the filled circles 
in Fig. 9 show the highest peak in each half-order. The fig- 
ure also shows published solar frequencies for I = 0-3 (Lazrek 
et al. 1997; Bertello et al. 2000; Chaphn et al. 2002). We see 
that the ridge of power follows the published frequencies very 
well, but again extends beyond them at high frequencies. It 
is interesting that our measurements at high frequency match 
up perfectly with the so-called pseudo-modes, also called HIPs 
(High-frequency Interference Peaks), which have been seen in 
smoothed solar power spectra up to 7.5 mHz (Garcia et al. 1998; 
Chaplin et al. 2001; Gelly et al. 2002). This may be relevant to 
the discussion of the physical nature of these HIPs: are they 
related to ordinary resonant p modes, despite having frequen- 
cies above the acoustic cut-off frequency in the atmosphere, 
and caused by small but non-zero reflectivity at the photosphere 
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(Balmforth & Gough 1990); or are they resonances between di- 
rect and reflected waves from a localized source (e.g., Kumar 
et al. 1990; Kumar & Lu 1991)? 

The points in Fig. 10 show the large separation as a func- 
tion of frequency for both a Cen B and the Sun. These val- 
ues were derived from Figs. 8 and 9, respectively, simply by 
doubling the differences between consecutive points along the 
ridges. The two curves are remarkably similar, except for the 
pronounced dip at 6mHz in a Cen B, which occurs in the re- 
gion of low S/N and which we suspect is not real. In fact, our 
tests show that this type of feature sometimes appears as an ar- 
tifact when the method is applied to short segments of the solar 
GOLF data (to match the observing window of a Cen B). The 
same tests also show that the method does measure Ai/ very ac- 
curately in the regions where oscillations have significant am- 
phtude. To illustrate this, the dashed lines in Fig. 10 show our 
best estimates, based on these tests, of the ±lcr uncertainties 
in Ai^. We conclude that this is a powerful technique for mea- 
suring the large separation over an extended frequency range. 
For a Cen B this will provide extra constraints on theoretical 
models while, for the Sun, we have been able to measure Ai/ 
to much higher frequencies than has been done previously. The 
S-shaped structure in the echelle diagram indicates a departure 
from the second-order fit used in equations (2)-(5), implying 
that cubic (or higher) terms are needed to describe the frequen- 
cies fully. As pointed out by Ulrich (1988), these higher-order 
coefficients might be useful in constraining the stellar mass. 

4. MODE LIFETIMES 

The scatter of the observed frequencies about the ridges al- 
lows us to estimate the mode lifetimes. In the case of a Cen A 
we did this by measuring the deviations of the measured fre- 
quencies from the fitted relations and comparing with simula- 
tions (Bedding et al. 2004). Here we adopted a slightly different 
approach that has the advantage of being independent of the fit, 
although it does still assume that we have correctly assigned n 
and / values. For each measured frequency j^„,/ (solid points in 
Fig. 6), except those at the ends of the ridges, we calculated the 
difference A„ ; between the measured frequency and that ex- 
pected from the positions of the two nearest neighbours on the 
same ridge, using linear interpolation: 

An- i'(„+An+),; + An+ ^(n-An-),; 
= An, + A„_ • 

Here, the two neighbouring modes lie in orders n + Am+ and 
n-Ati-, respectively. In most cases, A«_ and An+ are 1, but 
in some cases one or both is 2 or even 3. We therefore need 
to convert A„ / to a quantity that we can easily compare with 
simulations. We choose this to be the rms scatter of that peak 
about its expected position, calculated as follows: 

2 (A„.,)^(A«- + A«^)^ 

(An- + An+)2 + (A«-)2 + (An+)2' ^ ^ 

Averaging a over many peaks gives an estimate of the mode 
lifetime, but we must keep in mind that the finite SNR also 
introduces a scatter to the peak positions. To calibrate these 
two contributions, we carried out a large number of simulations 
(175 000), each with a single input frequency and each sam- 
pled with our observational window function (with the noise- 
optimized weights). We used the method described by Stello 
et al. (2004) to generate the time series of an oscillation that was 
re-excited continuously with random kicks and damped on a 
timescale that was an adjustable parameter (the mode lifetime). 



The other adjustable input parameter was the oscillation ampli- 
tude, while the mode frequency and the noise level were fixed. 
We made 100 simulations for each set of input parameters and 
from the resulting power spectra we measured the rms scatter 
in frequency of the highest peak and its mean SNR (discarding 
peaks with frequencies more than 4a from actual value). We 
repeated this for various values of the input parameters and the 
results are summarised in Fig. 1 1 . Each soUd line in this fig- 
ure shows the observed frequency scatter versus the observed 
SNR for a given value of the mode lifetime. As expected, the 
frequency scatter increases both with decreasing mode lifetime 
and with decreasing SNR. We stress that this figure applies 
specifically to our observing window for a Cen B and should 
be recalculated for other observing windows. 

The crosses in Fig. 1 1 show our results for a Cen B in two 
frequency ranges, centred at 3.6 mHz (the mean of all peaks 
below 4.0 mHz) and at 4.6 mHz (the mean of peaks above 
4.0 mHz). From these crosses we can determine the mode life- 
times. In order to check our results, we have also analysed seg- 
ments of the same 805-day GOLF series of full-disk velocity 
observations of the Sun used in Sec. 3.1. We divided this series 
into 100 segments, and imposed on each the a Cen B window 
function (with weights). The squares in Fig. 1 1 show the fre- 
quency scatter for each / = 1 mode as a function of the mean 
SNR (both measured over the 100 segments). We analysed 
twelve 1 = 1 modes (n = 16-27) and ignored the pairs with / = 
and 2 because they are less well resolved and interact via daily 
sidelobes. We see the well-known result that mode lifetimes in 
the Sun vary with frequency (the shortest lifetimes occur at the 
highest frequencies). Those with intermediate frequency have 
the best SNR because they have the highest amplitudes. 

For each measured point in Fig. 1 1 (two for a Cen B and 
twelve for the Sun), we used our simulations to infer the mode 
lifetime. The results are shown as a function of frequency in 
Fig. 12 using the same symbols, where we have expressed fre- 
quency in units of I'max, the frequency of maximum oscillation 
power (see Table 3 for values of I'max)- The diamonds in this 
figure are published measurements of the solar mode lifetimes 
(Chaplin et al. 1997), and we see good agreement with our val- 
ues. We also see that the typical mode Ufetimes for a Cen B, 
when considered as a function of v/vm-dx, are similar to those in 
the Sun. Our estimates are 3.3;l:o 9 d at 3.6 mHz and 1.9;I^4 ^ at 
4.6 mHz. 

4. 1 . Revised mode lifetimes for a Cen A 

The grey crosses in Fig. 12 show mode lifetimes for a Cen A, 
based on the observations analysed by Bedding et al. (2004). 
Here, we have re-measured the mode lifetimes using the 
method described above, which involved making a whole new 
series of simulations with the a Cen A window function in or- 
der to convert frequency scatters into lifetimes. The inferred 
lifetimes are 2.3^^;^ d at 2.1 mHz and 2.1-^% d at 2.6 mHz. 

We can compare these revised lifetimes for a Cen A with 
those we reported previously from the same data (see Table 3 
in Bedding et al. 2004): 1.4!*];|d at 2.1 mHz and 1.3![J:|d at 
2.6 mHz. The revised values are higher, although the Icr error 
bars do overlap. The reason for the change is that our previous 
calculation underestimated the contribution of SNR to the fre- 
quency scatter. We did of course include the effects of SNR, but 
did so by treating it as being independent of the scatter intro- 
duced by finite mode lifetimes. In fact, as our new simulations 
show, the two contributions are not independent. If they were. 
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the curves in Fig. 1 1 would not rise towards low SNR as steeply 
as they do. 

The important conclusion, as we can see from Fig. 12, is 
that mode lifetimes in a Cen A are not substantially lower than 
those in the Sun, although the value for the lower frequency 
range is still about 1.5 a below solar. 

5. AMPLITUDES AND NOISE LEVELS 

It is important to measure oscillation amplitudes in solar-like 
stars and to compare these with theoretical calculations (e.g., 
Houdek et al. 1999). It is also interesting to measure the back- 
ground noise from stellar convection, although in velocity this 
requires extremely precise measurements because the signa- 
ture is weak. For both these measurements, we have chosen 
to smooth the power spectrum heavily, so as to produce a sin- 
gle hump of excess power that is insensitive to the fact that the 
oscillation spectrum has discrete peaks. It is also useful to con- 
vert to power density, which is independent of the observing 
window and therefore allows us to compare noise levels. This 
is done by multiplying the power by the effective length of the 
observing run, which we calculated as the reciprocal of the area 
under the spectral window (in power). 

In Fig. 13 we show smoothed power density spectra for the 
Sun and four other stars: a Cen A and B, (5 Pav and /? Hyi. For 
the four stars, we used the most precise observations available: 
UVES observations for a Cen B (this paper), 5 Pav (this paper) 
and a Cen A (Butler et al. 2004), and UCLES observations for 
/? Hyi (Bedding et al. 2001). In all cases we used the raw veloc- 
ity measurements, before removal of any jumps or slow trends, 
since we are interested in measuring the total noise level. 

For the Sun, we used data from BiSON (Birmingham Solar 
Oscillations Network) and GOLF. The BiSON data comprised 
a 7-day time series with 40-second sampling from the Las Cam- 
panas station in Chile, kindly provided by W. Chaplin (private 
conmi.). The GOLF data comprised a 20-day time series with 
20-second sampling, kindly provided by P. Boumier (private 
comm.). Note that these have a higher Nyquist frequency than 
the pubhcly available GOLF data, which we used in Sees. 3.1 
and 4 but which are only sampled at 80 s. 

The dotted lines in Fig. 13 are fits to the noise backgrounds, 
based on the Harvey (1984) model of solar granulation. The 
Harvey model gives a convenient functional form, even in stars 
where the low-frequency noise has a strong additional contri- 
bution from instrumental noise. We discuss the noise levels in 
more detail below. First, we show that these smoothed power 
density spectra provide a powerful way to measure oscillation 
amplitudes in a way that is independent of mode lifetime. 

To do this, we first subtracted the background noise (dotted 
lines in Fig. 13) from each observed power density spectrum. 
We only included those parts of the spectrum that were at least 
twice the noise level. In order to calculate the amplitude per 
oscillation mode, we should then multiply by At//4 (where Az/ 
is the large frequency separation of the star) and take the square 
root. The rationale for this is that there are four modes in each 
segment of length Az/ (with / = 0, 1,2 and 3). However, we 
must keep in mind that modes with different angular degrees 
have different visibilities in full-disk observations, due to vary- 
ing amounts of cancellation. Based on the results presented by 
Bedding et al. (1996), we calculated the effective number of 
modes per Az/, normalised to the mean of the / = and l=\ 
amplitudes, to be 3.0. We therefore used this factor, rather than 
4, in our calculation. 



For IS.V we used the following values: 135 /kHz for the Sun; 
106 ^Hz for oi Cen A; 162 /zHz for a. Cen B; 56 /xHz for /3 Hyi 
and 93 /zHz for 5 Pav. The last of these is not a measure- 
ment, since none is available, and is instead derived from the 
following adopted parameters: M = 0.9Mq, L = I.SLq and 
Jeff = 5540 K. 

Our amplitude estimates are shown in Fig. 14, and the height, 
frequency and FWHM of the envelopes are given in Table 3. 
We see a number of interesting things from Figs. 13 and 14. 
Looking first at the Sun, the results illustrate very nicely that 
the solar oscillation amplitude depends on the spectral fine that 
is being measured (see Baudin et al. (2005) for a recent study of 
this phenomenon). The sodium Une used by GOLF is formed 
higher in the solar atmosphere than the potassium line used 
by BiSON, which is why GOLF measures higher oscillation 
amplitudes (Isaak et al. 1989). The actual height difference is 
difficult to estimate; Palle et al. (1992) quoted ~200km while 
Baudin et al. (2005) adopted ^60 km. Velocity measurements 
of other stars are made using a wide wavelength range so as 
to include many spectral lines, and the Doppler signal is domi- 
nated by neutral iron lines. Since these Unes are formed about 
400 km below the sodium D Unes (Elbe et al. 2001; Meunier 
& Kosovichev 2003), we would expect solar amplitudes mea- 
sured using the iodine technique to be less than those from both 
BiSON and GOLF. 

Unfortunately, there do not appear to be any published esti- 
mates of the solar oscillation amplitude using the stellar tech- 
nique. Here we present some previously unpublished observa- 
tions of the solar spectrum made using iodine referencing. One 
of us (J.T.W.) made observations of the Moon using the 0.6 m 
Coude Auxilliary Telescope (CAT) at Lick Observatory, which 
fed the Hamilton Spectrometer, a high-resolution (R = 60000) 
echelle (Vogt 1987). The CAT tracked a fixed point of uni- 
form surface brightness (Archimedes crater) for five consecu- 
tive nights near full moon, thus measuring the disk-integrated 
solar spectrum at night (a technique inspired by McMillan et al. 
1993). These velocity measurements allowed us to estimate 
the solar oscillation amplitude, which we include in Fig. 14 
and Table 3. The results support the conclusion that Fe I mea- 
surements give lower amphtudes than both GOLF and BiSON. 
They also place the solar amplitude between those of a Cen A 
and B, as would be expected given their stellar parameters. It 
would clearly be valuable to obtain more measurements of the 
Sun, in order to better calibrate the relationship between stel- 
lar and solar amplitudes. The method described here allows us 
to estimate amplitudes independently of mode lifetime and ob- 
serving window. Comparing the amphtudes of different stars 
with theoretical models is the subject of a future paper. 

We turn now to the noise levels in the various observations, 
looking first at low frequencies. The rise in power towards low 
frequencies seen in Fig. 13 is due to a combination of instru- 
mental drift and stellar background noise. Of course, it is very 
difficult to distinguish between these two, although in the Sun 
it is established the solar background is dominant. It therefore 
seems likely that the low-frequency power from a Cen A is also 
mostly stellar, given the power density is similar to the Sun 
(see also Kjeldsen et al. 1999). The same may also be true for 
a Cen B. We can certainly say that at 1 mHz, the granulation 
noise in both stars is no greater than is observed in the Sun by 
GOLF and BiSON. 

At the highest frequencies the noise levels are dominated by 
white noise from photon statistics. We can see that the power 
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density (which indicates noise per unit observing time) is low- 
est for BiSON and GOLF, followed by a Cen B (UVES) and 
a Cen A (UVES). In the last two columns of Table 3 we pro- 
vide an update to Table 1 of Butler et al. (2004), showing the 
noise per minute of observing at frequencies just above the p- 
mode envelope (2v„i^^ and also, where the sampling allows, 
at very high frequencies (llmHz). We should note that the 
power at high frequencies in the Sun (4-6 mHz), particularly 
in the sodium line used by GOLF, is dominated by solar noise 
that presumably arises from chromospheric effects, with only a 
small fraction being due to coherent p-modes. 

6. CONCLUSIONS 

Our observations of a Cen B from two sites have allowed us 
to identify 37 oscillation modes with / = 0-3. Fitting to these 
modes gave the large and small frequency separations as a func- 
tion of frequency. We also introduced a new method, involving 
smoothing in the \/S.v echelle diagram, that allowed us to trace 
the ridges of power, and hence measure the large separation, 
well beyond the central region. 

We inferred the mode lifetimes in two frequency ranges by 
measuring the scatter of the oscillation frequencies about a 
smooth trend, based on a calibration involving extensive sim- 
ulations. We found mode lifetimes in a Cen B, when consid- 
ered as a function of frequency relative to the maximum power. 



that are consistent with those seen in the Sun. We applied the 
same analysis to our observations of a Cen A and deduced re- 
vised mode lifetimes that are slightly higher than we previously 
published (Bedding et al. 2004). 

A limited set of observations of the star 5 Pav showed oscil- 
lations centred at 2.3 mHz with peak velocity amplitudes close 
to solar. Further observations are needed to determine the large 
separation and individual mode frequencies in this star. 

Finally, we also introduced a new method of measuring oscil- 
lation amplitudes from heavily-smoothed power density spec- 
tra. We estimated the ampUtude per mode of a Cen A and B, 
/3 Hyi, 5 Pav and the Sun, and pointed out that the results may 
depend on which spectral lines are used for the velocity mea- 
surements. 
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Fig. 1 . — Time series of velocity measurements of a Cen B, both before and after removal of slow trends. 
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Fig. 2. — Time series of velocity measurements of 5 Pav from UVES (diamonds) and UCLES (squares). The lower panel shows a 1.1-hr segment of UVES 
velocities, in which the 7-min periodicity is clearly visible. 
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Fig. 4. — Power spectrum of a Cen B from the combined UVES and UCLES data. Top panel: the noise-optimized power spectrum, in which weights were based 
on the measurement uncertainties, which minimizes the noise. Middle panel: same as top panel, but expanded to show only the central region. Bottom panel: the 
sidelobe-optimized power spectrum, in which weights were adjusted to minimize the aliases. The spectral windows are shown as insets and vertical dotted lines 
show our final oscillation frequencies, as listed in Table 1. 
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to the final frequencies, given by equations (2)-(5). 
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Fig. 6. — Echelle diagram of the final oscillation frequencies in a Cen B (filled diamonds), as listed in Table 1. Dotted curves show the fits to the frequencies, 
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Fig. 10. — The large separation as a function of frequency in a Cen B and the Sun, as derived from Figs. 8 and 9. The dashed lines indicate the ±lcr uncertainties 
on Au (see text). The solid line for a Cen B is the same as in the bottom panel of Fig. 7 and shows Az/q as calculated from equation (2). 
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Fig. 11. — Calibration of mode lifetimes for the a Cen B observing window, using the noise-optimized weights. SoUd lines are the results of simulations and 
show frequency scatter versus SNR for various mode lifetimes. Crosses show actual results for a Cen B in two frequency ranges, while open squares show results 
from GOLF data for twelve consecutive 1=1 modes in the Sun. 
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Fig. 12. — Mode lifetimes versus normalized frequency, where i^max is the frequency of maximum oscillation power. Values for a Cen B (black crosses) and the 
Sun observed by GOLF (open squares) are calculated from Fig. 1 1 . Two values for a Cen A (grey crosses) are calculated from a similar calibration of the results 
presented by Bedding et al. (2004). Open diamonds are published measurements of the solar mode lifetimes (Chaplin et al. 1997). 
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Fig. 13. — Smoothed power density spectra from velocity observations of the Sun and four other stars. The dotted lines are fits to the noise background. 
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Fig. 14. — Amplitude per mode for solar-like oscillations. These curves were calculated from those in Fig. 13 by subtracting the background noise, multiplying by 
A1//3.O and taking the square root (see text). We also show oscillation amplitudes in the Sun, measured from iodine-referenced observations at Lick Observatory. 
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Table 1 

Oscillation frequencies for a Cen B (/xHz) 



n 


1 = 


1=1 


1 = 2 


/ = 3 


1 n 
1 / 




iVjy. 1 ± U.y 






1 o 

Is 




J/z4./ ± U.9 






ly 


jjUo.o ± u.y 


Jj5l.9± I.I 


'3/1 c/: ^ _i_ 1 1 
3430.0 ± I.I 


ijlQ.i ±1.1 


20 


3466.9 ±1.0 


3544.9 ±0.8 




3685.6±1.1 


21 


3628.2 ± 1.0 




3778.8 ±1.1 


3849.3 ±1.3 


22 


3789.2 ±1.2 


3865.9 ±1.1 




4008.5 ±1.5 


23 


3951.1 ±1.2 


4025.9 ±1.2 


4102.0 ±1.2 




24 


4109.5±1.1 


4188.0±1.1 


4262.0 ±1.5 


4333.3 ±1.4 


25 


4275.7 ±1.5 


4351.0±1.4 


4425.4 ±1.5 




26 






4585.6±1.8 




27 


4598.4 ±1.6 


4670.3 ± 1.7 


4750.8 ±1.8 




28 




4835.4 ±2.0 


4912.4 ±2.0 




29 




4998.8 ±1.9 






30 


5085.7 ±2.2 


5155.6±2.1 






31 


5248.8 ±2.4 








32 


5411.4±1.9 


5489.0 ±2.3 







Table 2 

Frequency separations at 4.0 mHz for a Cen B (/xHz) 





161.50±0.11 


(0.07%) 




161.27 ±0.09 


(0.06%) 




161.48±0.17 


(0.11%) 




161.53±0.33 


(0.20%) 


Az/ 


161.38 ±0.06 


(0.04%) 




4.52 ±0.51 


(11%) 




10. 14 ±0.62 


(6%) 


Sfn 


16.73 ±0.65 


(4%) 


Do 


1.771 ±0.047 


(3%) 


£ 


1.477 ±0.011 


(0.7%) 



Table 3 

Amplitudes and noise levels for solar-like oscillations 







Peak amplitude 




FWHM 


Noise per min (ms ') 


Star 


Spectrograph 


per mode (ms~^) 


(mHz) 


(mHz) 


at 2 l^max 


at 11 mHz 


aCenB 


UVES 


0.085 ±0.004 


4.09 ±0.17 


1.98±0.14 


0.35 


0.32 


a Cen A 


UVES 


0.263 ±0.008 


2.41±0.13 


1.34 ±0.04 


0.49 


0.45 


/3Hyi 


UCLES 


0.432±0.016 


1.02 ±0.05 


0.54 ±0.05 


2.47 




5Pav 


UVES 


0.236 ±0.008 


2.33 ±0.09 


1.24 ±0.06 


0.82 




Sun 


GOLF 


0.308 ±0.005 


3.21 ±0.11 


1.76 ±0.02 


0.52 


0.20 


Sun 


BiSON 


0.233 ±0.006 


3.17±0.13 


1.62 ±0.04 


0.20 


0.12 


Sun 


Lick 


0.208 ±0.029 


2.95 ±0.31 









